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ABSTRACT 

This final report describes in detail the underlying theory and design 
of test instrumentation and computer techniques used to determine the 
power spectral density (PSD) of the Draper Laboratory Third Generation 
Gyro. Acceptance test results on the completed instrument (instrument 
build funded in part under this program) are presented as well as com- 
posite PSD plots determined from data taken at both Draper Laboratory 
and at Martin Company, Denver. The measured RMS noise over the 
frequency band 0. 0025 - 10 Hz was 0. 012 arc seconds (at Draper 
Laboratory). Results obtained at Martin Company indicate an RMS 
noise of 0. 009 arc seconds (0. 05 - 10 Hz). The results obtained showed 
the TGG performance was close to the proposed LST requirements but 
that additional development work would be needed to reduce the TGG 
noise to within more desirable limits. 
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CHAPTER 1 


INTRODUCTION 

The Large Space Telescope/ Control Moment Gyro (LST/CMG) 
test program attempted to establish the capability of low noise inertial 
gyros in precision stabilization and control of space vehicles. The 
basic approach used was two -fold: 

1) Establish gyro noise characteristics (18 IRIG Mod D and the 
Third Generation Gyro (TGG) ) in the 0.003 to 10Hz band. 

2) Determine gyro performance when used as an integral part 
of a simulated LST fine pointing system. 

The individual tasks included the setup and test of two 18 IRIG 
Mod D gyros in a back-to-back configuration, build and test (basic 
acceptance) a TGG from parts furnished GFE from the Air Force (SAMSO), 
establish TGG noise levels (at Draper Lab), and perform TGG noise measure 
ments and LST system tests at the Martin Company in Denver, Colorado. 
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CHAPTER 2 

DESCRIPTION OF TEST PROGRAM 

1. 18 IRIG Back -to -Back Testing 

The 18 IRIG Mod D gyro is a moderate performance rate integrating 
gyro designed and built at Draper Laboratory. It is a floated, viscous 
damped instrument with a gas bearing wheel, magnetic suspension, 
bridge type signal generator and soft iron torque generator. 

Two 18 IRIG's were mounted in a common test fixture (input axes 
parallel) and instrumented as shown in Fig. 1 . A brief description of the 
system shown in Fig. 1 is presented in the following paragraphs. A more 
detailed analysis of the electronics design and of the power spectral den- 
sity (PSD) techniques used are presented in Chapters 4 and 5. 

The outputs from the two 18 IRIG’s were amplified arid coupled 
through split winding transformers. The difference signal shown is the 
algebraic difference of the output signals from the two gyros. This signal 
was filtered, demodulated, and processed using PSD techniques. 

The difference signal eliminated any common mode output exhibited 
by the two gyros. The prime common mode signal was expected to be 
ground motion inputs to the two gyros. 

Also available for measurement and analysis were the individual 
gyro output signals. These signals were processed and compared to the 
difference output. The results obtained are discussed in Chapter 6. 

2. TGG Tests 

The Draper Lab TGG is a high performance, floated, viscous damped, 
rate integrating gyro. The instrument has extremely low drift **nd low noise 
and was considered an excellent candidate for the LST/CMG program. The 
TGG employs a gas bearing wheel, bridge type signal generator, and a soft 
iron torque generator. 
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Fig. 1 LST/CMG Teat, Gyro Noise Measurement System 




The basic acceptance test results for the TGG are presented in 
Chapter 3. After passing the basic acceptance test, the TGG was 
integrated with a system configuration similar to that employed with 
the 18 IRIG's. No back-to-back testing (TGG with 18 IRIG) was per- 
formed due to the high noise measured with the 18 IRIG (see Chapter 6). 

3. TGG Noise Measurements 

Noise measurements were obtained for the TGG at Draper Labor- 
atory and at the Martin Company, Denver, Colorado. Unfortunately, 
numerous computer related hardware and system software problems 
limited the overall success of the program and a complete characterization 
of the TGG noise could not be obtained during the initial contract period. 

Draper Laboratory requested a contract extension to allow Draper 
engineers the use of the contract deliverables on a separate program. 
During the extension period additional data was obtained at Draper Labor- 
atory which allowed a complete spectrum to be obtained for the TGG. 

This data is presented in Fig. 23 of Chapter 6. The RMS noise obtained 
in the band 0. 0025 - 10 Hz was 0. 012 arc second. 

Per NASA direction, the Martin test effort concentrated on obtaining 
a comprehensive TGG noise characterization. Due to lack of adequate 
funding the testing of the TGG on the Martin fine pointing simulator was 
not accomplished. 
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CHAPTER 3 


TGG ACCEPTANCE TEST SUMMARY 

Standard acceptance testa on TGG-GlA-221x were completed on 
4/12/74 prior to delivery of the instrument to the LST/CMG Program. 
Test results obtained are tabulated in Table 1 below. In addition to 
those results listed in Table 1. the following tests were performed and 
found to be within specified limits: thermal cycling (320° - 155°F), 
insulation resistance, signal generator and torque generator phasing, 
magnetic suspension tuning, float freedom, axial and radial centering, 
unbalance coefficient. 
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TABLE 1 

TGG ACCEPTANCE TEST RESULTS 


Float stops 

Signal generator null voltage 
Flotation temperature 
Gyro gain 

Torque generator sensitivity 
(per winding) 

Signal generator sideband noise 
(due to wheel unbalance) 

Total reaction torque 

Wheel excitation voltage 

Wheel operating power 

Average wheel power sensitivity 

Average wheel voltage sensitivity 

Wheel run-up time 

W r heel run -down time 


Numerical Value 

Units 

+4.0 - 4.7 

mrad 

3.0 

mV 

134.7 

°F 

130 

mV /mrad 

0.050 

dyn-cm/mA 2 

0.11 

mV (@533 Hz) 

<0.22 

dyn-cm 

45 

V 

10.8 

W 

0.015 

deg/h/W 

0.003 

deg/h/V 

30 

seconds 

110 

seconds 
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CHAPTER 4 


GYRO NOISE MEASUREMENT ELECTRONICS 


1. General Description (Fig. 2 ). 

The SG output of each gyro is amplified with a low noise preamp. 

By coupling the preamp outputs through split winding transformers, 
we obtain individual gyro signals and also a difference signal. 

A band pass filter provides selective amplification at the carrier 
frequency, and then a phase sensitive demodulator converts the signal 
to dc. A low pass filter provides the output for digitization by an A -D 
converter. 

The TGG and the 18 Mod D have, in effect, two torque generators. 

One TG, driven from a precision constant current source, is used to 
cancel the major portion of earth rotation and gyro drift. The other TG 
is used in one of two modes; 

(1) analog rebalance loop, or 

(2) constant dc current. 

When operating in mode (1), a track/hold circuit steps the output 
of a D-A converter to follow the analog torque current. When the track- 
ing circuit has settled, mode (2) provides a dc current at the proper value 
to prevent the float from drifting away from null . This technique allows 
automatic compensation for gyro drift. 
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Fig. 2 Gyro Noise Measurement System 




The following gyro parameters will be used below [refs. 1, 2J . 



TGG 

18 mod D 

units 

gyro gain 

1.2 

0. 516 


SG sensitivity 

15 

6.6 

V(RMS)/rad 

SG resistance at 6400 Hz 

66 

59 

ft 

float time constant 

750 

416 

jjisec 


2, Preamp (figure 3 ) . 

The preamp uses an LM381A integrated circuit op-amp refs. 3, 4 ] . 
This amplifier, marketed for high fidelity applications, has a noise voltage 
of 5 nV/ and a noise current of 0. 3 pA/ V5z at the carrier frequency 
(6400 Hz). 

The signal to noise ratio is maximum for a source impedance of 
16700ft, and this requires a 1:16 turns ratio transformer to couple the SG 
to the preamp. The most readily available transformers having this turns 
ratio also have so low a primary impedance that the SG would be excessively 
loaded. It was decided to postpone the purchase of a suitable transformer 
until its necessity is established. Failure to match impedances costs a loss 
of signal to noise ratio by a factor of 12. 

Because of the differential configuration, the amplifier noise is 
multiplied by 2. Using the gyro gain and the untuned SG sensitivity, we 
find the angle equivalent amplifier noise to be 0. 39 nr ad/ VS for the DSG 
and 2. 1 nrad/Vifr for the 18 mod D. For a 10 Hz bandwidth, the RMS noise 
would be 1. 25 nrad and 6. 6 nrad respectively. 

The preamp gain is 200, a factor of 100 coming from the closed 
loop gain of the op-amp, and a factor of 2 coming from the output trans- 
former. 

The grounding scheme for the SG and preamp is shown in figure 4. 
The input is kept balanced for maximum common mode rejection. 


9 


Note: Capacitors must be 
mylar or mica. 



200 SL 


II 


To 

bond posx 
4* i lie 







3. Band Pass Filter (figure 5). 

This filter is of the "negative immitance converter" type [ref. 5J . 
The transfer function is 


G(s") = - 


0<>./q) s 

s a + C'Jo/Q ) s ♦ < 


<0 


where G Q * 20, 

C0 0 = j/r X £>¥00 Hz , and ^ = C^o +/ )/2 * /0*S" . 

This function is plotted in figures 6 and 7 . 


Let us derive the effect of the filter on the demodulated signal. 


Signal: «- 

Carrier: C-OS Co 0 -£ - ^ ( Q*’*** 9 ^ + ^ 

Modulate: ^ tu, <> H ^ 

Filter: ^ [ G(S+i«>.) e (S+i "‘* + S(S-iM.) e**"^ ] 

Demodulate and normalize: ^fG(S+tw 0 )^G(S»ctOo)3 s si 


Hence the effect on the float angle is given by the transfer function 

HCs) « + G ( ^q)J (2) 

[ + 2Q (s/n.) +■ (sAo,,)* * Q Cs/oi.) 3 ^ 

I + t/Qls/v,) +(«Q%i) (S/w.) 1 +2Q (s/«.) 3 +Q i (s/^J * 

This function is plotted figure 8. 


4. Demodulator (figure 9 ). 


The demodulator uses an AH0014 integrated circuit analog switch 
[refs. 3, 4] driven by a digitally phase shiftable reference square wave. 
The LM308A amplifier was selected for its low drift (5 pV /°C), 
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Fig. 6 Gain of Band Pats Filter 
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Fig. 9 Demodulator 
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The demodulator gain is considered in steps. First convert from 
RMS to zero- to- peak, giving a gain of V 2 . Going from single ended to dif- 
ferential through the input transformer gives a gain of. 2. The amplifier 
has a gain of 1, but going from differential back to single ended through the 
amplifier gives a gain of 2. Also, there is a gain of 2/* in converting 
from zero- to- peak to average DC. 


A ripple filter is included, and the overall transfer function is 


(for RMS input and DC output) 


S-fU). 



where u> x = 2.*«r X / 00 Hz . 


5. Low Pass Filter (figure 10). 

This two- pole filter is used to reduce aliasing. Its gain and cut 
off frequency can easily be varied to suit test conditions. The transfer 

function is 


where uj t = 2w*/0 Hz, ani G* I A. -Oiz 
DSG electronics and G = 10 for the IS mod D electronics. 


(S+Wi ) 1 


[s) 


6. Power Amplifier (figure 11). 

This amplifier has a transistor output stage to allow higher loads 
than an integrated circuit op-amp could tolerate. The gain and output 
resistor are chosen to give a good closed loop transient response and a 
good operating range between saturation points. 

7. Track/ Hold Circuit (figure 1S». 

This circuit contains a 16 bit digital to analog converter (DAC) 
driven by an up-down counter. Although 16 bits is more than we need, 
the concomitant stability of its output motivates the choice of this DAC. 
It has a current source output in the range of ±1 mA capable of driving 
up to 1000 O. A 3PDT switch selects either track mode or hold mode. 
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Fig. 11 Power Amplifier 
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Fig. 12 Track/Hold Circuit 
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7.1 Track Mode. 

The TG is placed into an analog rebalance loop. An inverting 
low pass current filter (figure 1$ averages the TG current with a 20 
second time constant. The difference between the average TG current 
and the DAC current is the error input to the tracking loop. 

The error current is converted to voltage. The comparitor 
senses the sign of the error and steers the clock pulses into the up or 
down inputs of the counter so as to drive the error towards zero. 

7.2 Hold Mode. 

Constant current from the DAC flows throught the TG. The 
clock must be turned off before switching into hold mode. 


8. TG Configurations. 


8.1 The 18 IRIG mod D Gyro. 

The TG configuration is shown in figure 14. The torque expressed 
as angular velocity is [ref. 2, supplemented with communications from the 
gyro designersj 


With <*M49 jxrad/sec, the vertical component of earth rotation 
at CSDL, and with I s2 * 0, we get I p * I fll = 10. 4 mA, which is supplied by 
a precision constant current source. DAC current I s 2 is then applied with 
sensitivity 4. 7 (prad/sec)/mA. In hold mode the available current varies 
over a range of 2 mA, so the range of the track/hold circuit is 9. 4 ^rad/sec. 


8.2 Third Generation Gyro 

The TG configuration is shown in figure 15. The torque expressed 
as angular velocity is [ref. lj 
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Fig. 15 Torque Winding Configuration for the TGG 
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With I pE 3 I SE * 0. 7 mA and W * 49 pxad/sec, we get I po * 

Ig 0 * 11.4 mA, which is supplied by a precision constant current source. 
Then I pE = Ig E is supplied by the DAC. Since the DAC current is 0 to 
1 mA, we get a track/ hold range of 0. 38 j^rad/sec, 

9. Preamp Gain Adjustment 

In the dual gyro configuration, the gains of the two preamps 
must be adjusted so that the difference channel output shows no response 
to motion. This adjustment requires that the preamp outputs be in phase. 

SG tuning capacitors were not used because resonant circuits have rapidly 
varying phase shifts. There is no objection to the use of a tuning capacitor 
in the single gyro configuration. 

Unfortunately, the LM381A preamp requires a capacitor for proper 
operation |ref. 4] , and this introduces a phase shift. We have achieved 
an uncancelled motion sensitivity of 10 percent and expect to do better. 

10. Measurement System Gain and Transfer Function. 

The transfer function of the measurement system can be computed 
by multiplying together its components: gyro transfer function, which 
includes the float time constant (sec. l)j preamp gain 200; band pass filter 
(eq. 3); demodulator (eq. 4); and low pass filter (eq. 5), The frequency 
dependence is plotted in figures 16 and 17. The DC gains are tabulated 
below. Since they are based upon nominal gyro parameters, it is necessary 
to calibrate the system to obtain the actual values. 

The output of the low pass filter is digitized with a 12 bit analog to 
digital converter (ADC) whose input range is +5 volts. 



TGG 

18 mod D 

units 

Electronic gain (RMS input. 




DC output) 

14400 

144000 


Measurement system gain 

260000 

490000 

V/rad 

ADC range 

+19 

+10.2 

f*rad 

ADC quantization 

9.4 

5.0 

nrad 

Note: Hased upon nominal gyro parameters. 
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Actual DC gains (based on nominal gyro parameters), 

TGG 260000 V/rad 

TGG electronics only 14400 

18 Mod D 490000 V/rad 

18 Mod D electronics only 14400 
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CHAPTER 5 

POWER SPECTRAL DENSITY TECHNIQUES 


This paper is an introduction to Fourier transforms and power 
spectral density techniques in the measurement of signals and noise. 
Section 1 provides the mathematical basis and includes explanations of 
aliasing, Fourier tranforms, power spectral density, and the fast 
Fourier transform. Section 2 describes the averaging which must be 
used when measuring the power spectral density of noise. The im- 
portance of power spectral density is explained at the end of this section. 
Finally, section 3 gives some helpful computer methods. 

1. Mathematical Basis. 

1.1. Positive and Negative Frequencies. 


A signal of a single frequency f is given by 


CL 


ZTift 

c 


•f 


* -lirift 

cu e. 


( 1 ) 


The two terms in (1) are at positive and negative frequencies. The 
* 

amplitudes g. and g must be complex conjugates of each other in order 
for x(t) to be real. Equation (1) is equivalent to 

id)- 2la.lcos(tirfe +#), m 

so that the use of positive and negative frequencies is just a convenient 
way to specify phase information. Mathematical frequencies f and -f 
always appear together in a signal of physical frequency f . 
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1, 2. Sampling and Aliasing 


( 3 ) 

be sampled at interval At. The sampled values are 

- C») 

It is seen that for any integer m, if the frequency f is replaced by f + m/At, 
the sampled values in (4) remain unchanged. In particular the two signals 

m / 

C^ a ") e 7irlft 

yield identical samples. 



Let the signal 


= £ 




In order to insure a unique relation between the continuous signal 
and the discrete samples, we must restrict the physical frequency to the 
range 

0C-T< 2 ^ * • (0 

The mathematical frequency is then restricted to the range ~* max <? < f m<ot * 
an interval of length 1/At. Usually * max is called the Nyquist frequency. 

If frequencies greater than f max are present, and the samples are 

processed anyway, then these higher frequencies will appear as a spurious 

signal below f in the manner indicated by equation (5). This distortion 
max 

is called aliasing. 


Example: Suppose At * 0. 01 sec, so that f max * 50 Hz. Suppose 
an f * 60 Hz signal is present. An integer multiple of 1/At « 100 Hz 
must be added to f to bring it into the range -50 <f <50. Thus 60 + (-1) (100) *-40, 
So a 40 Hz spurious signal will appear. Furthermore the complex amplitude 
at+60Hz is added to the already existing complex amplitude at -40 Hz. 
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1. 3. Discrete Fourier Transform 


In what follows, use will be made of the relation 

j«o ( 0, if mod » . 

Equation (7) is easily proved by representing the summands as unit vectors 
in the complex plane. 


Let x^, k = 0, 1, .... n-1 be sampled values and define the Fourier 
trainsform by 

n-l , :w 

u _ i 


i/j * I e 
J fc-0 


It is seen that 


«J»i> 


so that yj is periodic in j with period n. It will sometimes be convenient 
to regard the range O^j ^n**l as the fundamental period, and sometimes 
the range 

J. + I C j C ^ > «P»r n evaa 

... } do) 

2 ~ < j $ ^ , <V n odd . J 


Por n add . 


The inverse Fourier transform is 

. ji 

I 4 2l»t J* . % 

" ?,4 e * «/. . (ll) 

J*0 d 

To verify this, substitute (8) into the right hand side of (11). Using (7), 
the right hand side simplifies to x^. Since (8) and (11) are linear transforms 
on an n- dimensional vector space whose product is the identity, they are 
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J I L 


Inverses of each other. 

Equation(ll) expresses the signal x(t) as a superposition of single 
frequency signals. If we wish to put equation (11) into the form 


2iri *P. k At 

J 


= ^(kAt) ^ J t (tz) 

we must have J 

f i k M * A + »*? (U) 

J h ^ 

for some integers m. (because the exponential function has period 2vi). 

J • 


r _ J 
+ j ~ r>4t 


*j 


k At 


Oi) 


Since r cannot depend upon k, we must have m^ = mk, with m an integer, 
and so I m. 

= + -2- . Os) 

J hAi At 


According to the restriction (6), we set m * 0, and in (12) sum j over the 
range given by (10). 


Using the notations 


Af~ 


l 

h At ’ 


a c j Af ) - uj » 


we put (12) in the form 


■*«) = -prZ a 0 a^) e 




Ut) 


07) 


(If) 


which is manifestly a superposition of single frequency signals. The 
summation range is given by (10), and is approximately -n/2 < j <n/2. 
Therefore the maximum frequency is (n/2) / (nAt) * * max » l* 1 agreement 
with (6). 


1.4. Degrees of Freedom . 


Reality of requires 
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and in part cvlar y^ is real. Periodicity of y^, equation (9), requires 
that y n j 2 he real. 


The sampled values are n real numbers. The Fourier coefficients 
y^ also have n degrees of freedom: yg, which is real; if n is odd, then 
(n-l)/2 complex numbers for 1 ^(n-l)/2; if n is even, then n/2-1 

complex numbers y^ for 1 ^n/2-1 and y^g* wh ich is real. 


It is easily seen why y^g has one rather than two degrees of freedom. 
At this frequency the sampling occurs at 180° phase intervals. If t * 0 
is one of the sample times, then the signal sin(2xft) will always be zero 
when it is sampled. 


1. 5, Power Spectral Density . 


The power in the signal x(t) is defined by 


P = T C dt *7^ i * (k Ai) * M Cio) 


l-i 


7T 2 * 


!*0 


Substitute for from (11) and use (7) to obtain 

5 l»ji‘ 


_ I 


Z ** • 


h* Afi \) 

For a continuous frequency domain, the power spectral 
density is defined by 


J - 1 •» a# 


»0 


( 22 ) 

(22) 


( 2*0 


S3 


J 


1 


1 


X 


If we now regard (23) be being summed over the range (10) and regard it 


for j / 0, n/2. Using (19) and (16), 


PSDtjaO- 


to (24). 

we have 


1 

-(IV- 

*!Vil l ) . 

(16). 



2 At 
n 

l*jl\ 

°> 2 

h 

h 1 

V 2 , 

j-o,i 


(2r) 


( 26 ) 


The cases j * 0 or n/2 are special since there is only one 
degree of freedom at frequencies f= 0 and f = f max compared to the usual 
two degrees of freedom at other frequencies. The case f = f max 
can be ignored because if the component at this frequency is not 
negligible, then there will usually be non- negligible components at 
slightly higher frequencies, and these will cause aliasing. The case of 
zero frequency may or may not be important in particular applications. 


1. 6. Filtered Signals . 


If the signal x(t) is the input to a linear system with transfer function 
G(s), and if u(t) is the output, then u(t) has a Fourier transform V(f). 

It is clear that the relations between the two Fourier transforms and the 
two power spectral densities are 


Vf-P) = G(2iri-T) a f-f), (») 


PS 0 tt «). I 


Example: An integrating gyroscope provides a direct readout of 
angle. By use of an analog torque to balance loop, we can get a readmit 

of angular velocity. The power spectral density of the instrument noise 

2 

is given in the former case as PSDg(f) in units of radian /Hz, and in 
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the latter case as PSD^ff) in units of (radian/sec) /Hz. These two repre^ 
sentations of the gyroscope noise are related by 

PSD 4 (f) = (2*f) 2 PSD fi (f). 
u o 

1. 7. RMS V alues . 

Prom (21), (22) and (25) 

% 2 = p = £ PSDCj A-f) M. m 

If x(t) is filtered with an ideal band pass filter which only passes frequencies 
between f, and f 0 , then it follows from (28) that the RMS value of the 
output is 2 

Cf.^) = I • PSbCjA-f) Af - (30) 

In the continuous case, 

2 f 

Vzhs C-f./O *• L 2 PSD(-C)d-f- »') 

J $\ 

By allowing fj and fg to become close together, we see the physical 
significance of the power spectral density: PSD(f) Af is the amount of power 
present in the signal in a bandwidth Af about frequency f. 

If x(t) is measured in units X (e. g. volts), then the power spectral 

2 , 

density is measured in units X /Hz. 

1. 8. Effects of Finite Frequency Resolution. 

We can think of a signal as possessing a power spectral density, which 
we want to measure. If the duration of the test is T* nAt, then the 
measurement has a frequency resolution Af * 1/T. 
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If the signal has a broad band power spectrum, then for sufficiently 

« 

small Af, a further decrease will not change the measured values of PSD(f) 
On the other hand, if the signal contains a finite amount of power at a 
single frequency f Q , then the true power spectral density must have an 
impulse function at fg because the integral in (31) takes a jump in value 
when one of the limits of integration passes through f Q . In this case, 
as Af decreases, the measured PSD(f) will increase indefinitely in a 
decreasing neighborhood of f q. 


As an example, let us consider a single frequency signal 

<«*)-. ci + a* e- 1 2 ^* 

Let j Q be the nearest integer to f Q /Af. Then we expect that the signi- 
ficant values of y.. will occur for j«±j 0 . Assume that 1 and 
n/2 - j 0 »l, which means that the positive and negative frequency parts 
are separated by a large multiple of Af. A straightforward derivation 
then yields. 


PSD CjAf) 


2 q 2 

A f 



( 33 ) 


If f ft is not an integer multiple of Af , then the measured power 
w . 2 

spectral density will be spread out with side lobes falling off as (j -jg) . 


1. 9. Hanning Window 

The (j - j 0 )~ Z spread is characteristic of a discontinuity in the 
data. From a mathematical point of view the data x^ are periodic 
with period n, and therefore there is an end point discontinuity between 

Vl and *n = *V 


The end point discontinuity can be eliminated by multiplying the 
data by a window function which smoothly falls to zero at the end points. 
The Hanning window is 


w* * VI ( 1 ~ 


Co-0 


2wk \ 
n ) 


(3V) 


the factor V 2/3 is needed to keep the average power in ^ 

the same as in the original 2^. 


It is usually easier to apply the Hanning window in the frequency 
domain. If z • and are the Fourier transforms of and respectively, 
then for an arbitrary window 


n-» 


»i * I ' 


4=o 

and for the Hanning window 


( 3 $) 



( 31 ) 


In (35) and (36) the periodicity of the Fourier transform is implicitly 
used. 


Considering again the signal in (32) and the same approximations 
used above, we get for the power spectral density with Hanning 


PS D(jAf) 


4 a* 



Now the side lobes fall off as (j-j 0 )* 6 . This faster fall off Is a 
consequence of the continuity of the signal and its first derivative at the 

end points. 


( 37 ) 


The price paid for smaller side lobes is a more spread out center 

lobe. 
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1.10 The Fast Fourier Transform (FFT). 

The FFT is an algorithm discovered by Cooley and Tukey 
(Mathematics of Computation, Vol. 19, 1965, p. 297), which 

enables the Fourier transform to be computed in a reasonable amount 
of time. It will be assumed that n is a power of 2. 

Let x^, k * 0, . . . , n-1 be the given sequence, whose Fourier 
transform y., j = 0, . . . , n-1 is sought. Let x^, k = 0, ... n/2-1 be the 
sequence of even terms, = x 2k , and let y^, j = 0, . . . , n/2-1 be its 

Fourier transform. Let x" k , k ■ 0 n/2-1 be the sequence of odd 

terms. x M k = x 2fc f r and let y”.. j = 0 n/2-1 be its Fourier trans- 

form. 


9 :- 


From the definition (8), 

fc/2 ” -2*; il2S 

e h 


^2k 


+ ^ e ‘X 2 k* 

k=0 


= 2. e »/* % k +• ^ 

k«’o 


* n/a-» . jV 

h / & ^ 


k-o 


/ -2*i i. 

» y j + e * y 


CBS 


When j ^n/2, y.’ and y"j are defined by their periodicity, equation (9), 
where here the period is n/2. Equation (38) can be written as 

y.= fj + e'^i 


u" . 
J J 


1 




To obtain y/ and y/'. apply equations (39) recursively, with n/2 

in place of n, y'^ or y'\ in place of y^ and appropriate subsequences in 

place of y’. and y" . Since n is a power of 2. this recursion can be applied 
J J* 

until n*l, where the Fourier transform of a single number is that same 
number. 

Similar algorithms can be derived for prime numbers other than 2. 
For each distinct prime factor of n, a separate set of equations is needed, 
(39) being the set for the prime number 2. However, the FFT algorithm is 
most efficient when n is a power of 2. 

Computing the transform directly from (8) as a matrix product 
2 

requires n multiplications. Using the FFT, equation (39), requires 
n loggn multiplications. If we take as a typical case n=2048 and assume 
1 msec to do a complex multiplication, then use of the FFT reduces the 
computation time from 70 minutes to 23 seconds. 


2. Analysis of Noise . 

The values of the signal at particular instants of time are random 
and satisfy some probability distribution. If the noise has non- zero mean, 
lump this mean with the non-random part of the signal. 

We assume the noise is stationary, i. e. , that the statistical properties 
are unchanged by a shift in time. 

2. 1. Gaussian Noise . 

The following analysis is rendered simple by assuming that the 
noise Is Gaussian. The values of the n samples Xq.Xj, . . . , x Ra>1 are 
given by a Gaussian joint probability distribution. Since the Fourier 
transform, equation (8), is linear, the real and imaginary parts of all 
the y^ also have a Gaussian joint probability distribution. 
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For a fixed frequency jAf examine the two-dimensional probability 


distribution of y,. R and y^, the real and imaginary parts of y^. It is 


Gaussian with covariance matrix. 

C = 


« K 

L P®*®! 


r «■„<?! 

<Ti* 


(«<£>) 


For a time shift corresponding to a phase shift 0, 

* 


i($ 

yj is multiplied by e r . This is a rotation of the y. R and y^ axes in 
the y^ plane. Matrix C is left unchanged by all such rotations of axes if 
and only if 

<:.■»* • 

0 


■r 


*j- 


(Hi) 


Thus, if the noise is stationary, the real and imaginary parts of y. are 
independent. 


The y. R , y^ probability distribution is 


^ V-— - — (3jI ' Ul) I) , ( 12 ) 

\ 2. cr 2 / 

where u R , Uj are the non- random mean values. With y.j held constant, 
the probability distribution in y^ R has mean u R and variance Therefore, 
using an overbar to denote mean value. 


V - 


- (5j 0 * + ** .« ««* + «■* . 


Similarly, 


y • * 


s U 


i * cr* 


Adding (MS) and (4M). 


yj |* * +«i l * ia-*. 


(n) 


(it) 


(‘is-) 
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Now, u R and Uj depend upon the relative phase between the signal 
and the sample times. 


U 


M Uj * U ^ ) 


hjr = u 2 + 2<r 2 . 


(Hi) 


.2 


The quantity |y^ | is a random variable with mean given by (46). 

To obtain its standard deviation, consider far simplicity the two limiting 
cases u » <r and u = 0. 

In the former case, use the techniques of error analysis. 

4 ( lyjl 2 ) * 2 y jR Ay jR +2 Jjjj A 3^ . C*P) 

Combine the errors as independent and use y^ R *<u R , y^ <*Uj, Ay^ R = Ay^ *<r. 

*0ty 2 ) _ = 2£_ ^ 

i 2 " 2 LX 

I */ • «• 

2 

If u * 0, then from (26), |y^.| has a chi-square distribution with 
2 degrees of freedom, and this distribution is 

0 




2 <6 


_ 2 


20 " 


l e J ? >0 • 


(m«0 


The mean and standard deviation are both equal to 2o^. 

Summary of this section: If, at a given frequency f, the sig! 
consists of a non- random component of amplitude u and a Gaussian ran<!5m 
component having standard deviation <Mn each of its two degrees of 
freedom, and if the power spectral density. PSD(f), is calculated by 
equation (26), then PSD(f) is a random variable for which the ratio 
(standard deviation) /(mean) has t))e value 1 for u « e, the value 
2c/u for u »<r, and varies smoothly in between. 
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Conclusion of this section: In order to obtain a reliable estimate 
of PSD(f), it is necessary to make repeated measurements and average. 

If N measurements are averaged, the standard deviation in the estimate 
will be reduced by a factor of If we want the ratio (standard deviation) 
/(mean) for the estimate to be € , then in the: most pessimistic case (u=0), 
we must have 



Note that this conclusion depends to some extent on the assumption of 
Gaussian noise. 

2. 2. Importance of Power Spectral Density in the Description of Noise in 
Linear System s . 

Power spectral density is an essential tool in describing noise in 
linear systems. If we know the PSD of the input noise, equation (28) 
gives the output noise PSD, and equation (31) gives the RMS value 
of the noise.' Alternatively, if we know the output noise PSD and the 
transfer function, then we can calculate the input noise PSD. When in- 
dependent noise sources are added, the PSDof the sum noise is 
conveniently the sum of the individual PSD ’s. 

The significance of the PSD is that it carries the necessay 
frequency dependent information. If we only knew the RMS input noise 
to a frequency depdendent linear system, we would be unable to compute 
the RMS output noise. 

Example: A resistance R and capacitance C are in parallel. What 
is the RMS thermal noise across the parallel combination? 
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The resistance generates thermal noise with power spectral density 

pSD', n (-f)= 


-23 . 


in units of volt 4 /Hz. Boltzmann constant: k s 1.38 x 10 joule/ K. 


The gain of the system is 


\ 




The output PSD is 


pso. ut (f ) = -iiLg— 

i i — • r i » 


VkTR 


The RMS output noise is 

51 r * ifkTR 
RMS 


/ / + airt-T-P / 1 1 + 


r co 
J o 


_ ^ s tfkTR _3 __ feT 

0 I + (*rvf)‘ ~ ‘iirr Q ' 


For T * 300 °K, C = l^F, this becomes V RMg - 65 nanovolts. 

3. Computational Techniques . 

The most important computational technique is the FFT, described 
in section l.fC. 

3. 1. Efficient Use of Memory. 

The FFT algorithm can be programmed in such a way that each 
time equation (3? ) is applied, the new values are stored back in the same 
memory locations occupied by the old values. To perform an FFT on 
a sequence of length n, one needs 2n words of memory: n for the data 
(real part of x^), n of zeroes (imaginary part of x^). One also needs . 
n/2 words to hold the power spectral density and n/4 words for a table 
of sines. 
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By separating equation (3?) into its real and imaginary parts, one 
can carry out the FPT algorithm using only n words of memory. Furthermore, 
this can also be programmed so that at each step the same memory locations 

are reused. 

Another computational trick allows two Fourier transforms to be 
taken in the same FFT operation. From (19), if are real, the real 
part of y. is an even function of 3 and the imaginary part of y^ is an odd 
function of 3 . Similarly if x k are imaginary, the real part of y^ is an 
odd function of j and the imaginary part of y^ is an even function of 3 . 

Therefore, two sets of data can be stored as the real and imaginary parts 
of Xj^, and their respective Fourier transforms can be separated at the 
end of the calculation. 


Suppose the two sets of data are independent sets of samples of the 
same signal. Applying equation (26) for j j 0, n/2, we see that the res! 


data gives 

psd,.. 

hk'+t. 

, (S') 

and the imaginary data gives 


. («) 

PSD- „ L [ a / 

V. 2. / J 


If these two measurements are averaged, we get, after some easy 


algebra, 

psD j- tH |a J |! ' 

la-il 1 )- 

C Si) 


This allows the power spectral density averaging technique to be done 
two at a time. 


!. Digital Filtering to Reduce Allaslu 

Suppose the analog to digital converter can convert at time interval 
AT, but we wr.nt a sample time 

At - r at . <"> 


Instead of letting the computer to idle for time At, we can average 
r data. Thus we program a digital filter whose input is = x (kAT) and 
whose output is 

Uj» u.(jAt) = X £ (ss) 

Kao 


If the input is at frequency f, 

' iiri f ( K AT) 

*-k * e 

then the output is 


The filter output covers a frequency range ^ ' 


(Sfe) 


J I 'n£at\ I 6 


<S 7 ) 


(St) 


and because of aliasing, this frequency f receives components from 
the set of frequencies 

f, -l f *irl> 

where m is an integer. If the components at these different frequencies 
are uncorrelated, then the power spectral densities are additive, but 
with transfer function weighting given by (57). 
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ttie-O* 




•LSSflfar «■**•»• 

In the limit r-* oo , AT-* 0, At held constant, we get 


(6C>) 


. * xi^. Z (ir-f At) 

p sb «w>« I ISZo- pst> »(K‘S-0- 


M a -00 


( 6 /) 


Although this simple digital filter does not do a terrific job of 
reducing aliasing, it is much better than none at all. A more complicated 
filter would slow down the computer and require AT to be larger. Alter- 
natively, a hardware digital filter could be constructed. 
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CHAPTER 6 


TEST RESULTS 


l. 18 IRTO Mod D Tests 

Two 18 IRIG Mod D gyros were set up in a back-to-back configuration 
as shown schematically in Fig. 1. A power spectral density (PSD) was deter- 
mined for each gyro, as well as for the difference output. These results 
are shown in Fig 18-20. 

The magnitudes of the PSD's determined were much larger than 
expected. The prime source of the noise appeared to be in the gyro wheel. 
Substantial noise also coupled from the magnetic suspension. 

The noise coupled from ground motion would be expected to cancel 
in the difference output. However this output was comparable in size to 
those obtained with the individual gyros, indicating that any noise contri- 
bution from the ground was masked by the random noise emanating from 
the gyros. 

Fig.21 shows the PSD obtained with the input to the gyro electronics 
shorted. The noise is well below that measured with the gyros and has an 
RMS value of approximately 1,5 milli-arcseconds. 

2. TGG Test Results at Draper Lab 

The TGG was tested at Draper Lab using an electronics design 
patterned after that used with the 18 IRIG. With the input to the electronics 
shorted the equivalent TGG angle noise is presented In Fig. 22. The noise 
level is somewhat lower than that obtained for the 18 IRIG. having an RMS 
value of approximately 1 milli-arcsecond over a frequency band of 0. 01-10hz. 
The difference is due primarily to the higher TGG SG sensitivity. Fig. 23 
presents a composite PSD obtained for the TGG over the band from approx - 
imately 0. 0025. to 10 Hz. The RMS noise over the total band is It milli- 
arcseconds, 

3, LST/CMG Tests at Martin Company 

LST/CMG tests at Martin Company were limited to TGG 
noise measurements on one of the Martin test piers. After the initial setup 
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Fig. 18 18 Irig Power Spectral Density, A Gyro, Normal Conditions 














of equipment the TGG was position IA vertical and stabilized with constant 
dc current applied to the gyro torquers. Power spectral density measurements 
were made over a frequency range from approximately 0. 05 Hz to 10 Hz 
(Fig. 24, 25). A composite PSD graph with averaging over discrete frequency 
bands is shown in Fig. 26 . The RMS noise over the band from 0. 05 - 10 Hz is 
approximately 8 milli-arcseconds. 

During the contract extension period, data taken on the test pier with IA 
horizontal was recorded on magnetic tape and later was reduced at Draper 
Laboratory. Fig. 27 presents a computed PSD over the frequency range from 
approximately 0.05 to 10 Hz. Some averaging was performed resulting in 
a more accurate, smoother plot. The data is in good agreement with that 
shown in Fig. 26 . 

The occasional spikes shown in the plot are believed due to gyro tem- 
perature controller noise. The RMS noise for the band shown is 0. 009 arc 
second. Without including the larger spikes the RMS noise from 0.05 to 
10 Hz is 0. 005 arc second. 

Fig. 28 presents a PSD obtained from the same data tape for the 
frequency range from ~ 0. 01 Hz to 0. 5 Hz. The results agree well with 
data taken at Draper Laboratory (Fig. 23 ). 
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CHAPTER 7 


CONCLUSIONS 

Based on Draper Laboratory data (Fig. 23 ) the TGG noise over a 
frequency band from 0. 0025 to 10 Hz is 12 milli-arcsecond. Additional 
data recorded with the TGG on the test pier at Martin Company and later 
processed at Draper Laboratory shows that the true gyro noise in a 
seismicly quiet environment may be somewhat less than this value: i. e. , 
0. 009 arc second over the 0. 0025 - 10 Hz band. 

The test results showed that the TGG is compatible with LST require- 
ments, though the frequency band over which the gyro would maintain 
control would need to be reduced if a maximum 0. 005 arc second RMS 
noise were required. 

To reduce the gyro noise to an RMS value near one milli-arcsecond 
will require changes in the instrument; i. e. shrouding the wheel and/or 
making changes in instrument gain and suspension mechanization to 
minimize output axis coupling from the wheel and other instrument error 
torques. 
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